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Abstract 

For the Davey-Stewartson I equation, which is an integrable equation in 1+2 dimen- 
sions, we have already found its Lax pair in 1+1 dimensional form by nonlinear constraints. 
This paper deals with the second nonlinearization of this 1+1 dimensional system to get 
three 1+0 dimensional Hamiltonian systems with a constraint of Neumann type. The 
full set of involutive conserved integrals is obtained and their functional independence is 
proved. Therefore, the Hamiltonian systems are completely integrable in Liouville sense. 
A periodic solution of the Davey-Stewartson I equation is obtained by solving these clas- 
sical Hamiltonian systems as an example. 



1 Introduction 



The Davey-Stewartson I (DSI) equation is a famous 1+2 dimensional integrable equation de- 
scribing the motion of water wave |J. It has been discussed by various methods. Soliton solu- 
tions were obtained by inverse scattering method ]7j, [13], Backlund transformation [I], binary 
Darboux transformation [12], nonlinearization method [O, [18] etc. Almost-periodic solutions 



were also obtained 
In 



0- 



17, 18], this 1+2 dimensional problem was nonlinearized to an essentially 1+1 dimen- 



sional linear system (|l|) where all the differentials are separated. This system is very useful to 
get localized soliton solutions by Darboux transformation in 1+1 dimensions. 
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On the other hand, many 1+1 dimensional integrable systems can be nonlinearized to 1+0 
dimensional (or called finite-dimensional) integrable systems ]2j, |], [15|, [16], and the idea of 
nonlinearization was proposed by Cao [Q. Some 1+1 dimensional problems were completely 
solved and the periodic or quasi-periodic solutions were obtained. 

In 0, the KP equation, which is 1+2 dimensional, was nonlinearized not only to 1+1 
dimensional j|, §, but also to 1+0 dimensional integrable Hamiltonian systems. 

In the present paper, we show that the 1+1 dimensional system obtained by nonlinearizing 
the Lax pair of the DSI equation can also be nonlinearized to three 1+0 dimensional Hamil- 
tonian systems. We find a full set of involutive conserved integrals and prove their functional 
independence. Therefore, these systems are completely integrable in Liouville sense. As an 
example, when the number of the eigenvalues is two, we solve the systems directly to get a 
periodic solution of the DSI equation. 

It is well known that the DSI equation has a Lax pair in 1+2 dimensions. In [17, 18], a new 
integrable system was presented, which is essentially 1+1 dimensional, since all the differentials 
are separated. This system can be written explicitly as 



$„ = V<t> 



iA u if 
— u — iA — ig 

\ i/ -i<? J 



V 



-2iA 2 + i|w| 2 + it>i 
2uX + iu y 

-2ifX + 2fy 



—2uX + iu y 
2iA 2 -i|w| 2 - 
2igX + 2g y 



1 iA if y 
iA ig $ 

-2ifX - 2f y 
iv 2 2igX - 2g y 



•2KI/P 



(1) 



\gn J 



Here u, f and g are complex functions, v\ and v<i are real functions. 

Its integrability conditions § xy = $ yt = $ ty and $ xt = $ te consist of the following 
three parts. 

(1) DSI equation 

-iu t = u xx + u yy + 2\u\ 2 u + 2(vi + v 2 )u 

V\, x ~ Vl,y = V 2 ,x + V 2 , y = -(I^Hx- 

(2) Standard Lax pair of DSI equation 



1 
-1 



u 
-u 





(i o N 




f u\ ( 


F t = 2i 




| F xx + 2i 


\F X + i 




V° ~\ 







\u\ 2 + 2vi 

— U X + Uy 



u x + u 



(3) 



y 
2v 2 



where 




(4) 
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(3) Nonlinear constraint 



FF* = - 



v\ u x 



(5) 



u , 



i' 2 



Therefore, the nonlinear equations we will consider consist of the DSI equation, its standard 
Lax pair in 1+2 dimensions and the nonlinear constraint. As soon as this problem is solved, 
we get the solution of the DSI equation, although only a part of the solutions can be obtained 
by this method. In |18| , localized (M, N) soliton solutions were obtained from the Darboux 
transformation for ([!]). 

As many other 1+1 dimensional problems, here we want to find the nonlinear constraint for 
([]]) and perform nonlinearization again to obtain 1+0 dimensional Hamiltonian systems. 



2 Nonlinearization and Hamiltonians 

In order to get a nonlinear constraint which is compatible with all the x, y and i-equations in 
(H), we must consider the ^/-equation first. 

Notice that if $ is a solution of ([I]) for real A, then 

= -V T § $ x = -U T § $i = -W T $>. (6) 

Suppose $ = (0i, 2 , 03) T , then we may choose (i0i, i0 2 , i03) T to be the corresponding conjugate 
coordinates. 

The first element of the Lenard sequence corresponding to w — {u, —u, if, if, —ig, —ig) is 

\U 

G = (0, 0, if, if, ig, ig). (7) 
By computing the variation of A 0, we have 

5X/5UJ = (i0l0 2 , 10201, 10103, 10301, 10203, 10302)- (8) 

Now take N distinct eigenvalues Xj e 1R (N > 2). Suppose the corresponding solution of 
the Lax pair for A = Xj is (0ij, 2 j, 03j) T - Then define 

A = diag(Ai,---,A 7V ) = (0 j i,---,0, JV ) T (j = 1,2,3). (9) 

Each <l>j is a column vector, and A is a constant real diagonal matrix. 

We choose (4>ij, 02j, 03j, i0ij, i02j, i03j) as the coordinates instead of the real ones Re(0jfc) 
and lm(0jfc) (1 < j < 3, 1 < k < N). The standard symplectic form of 1R 6N is given by 

lu = 2 dlm(0j a ) A dRe(0j a ) = ^ i d0j a A d0j a . (10) 

l<j<N i<i<v 

l<o<3 l<a<3 

Denote (Vi, V2) = V*V2 for two column vectors V\ and V^. 
Define the nonlinear constraint 

G = E 5 V^ (11) 

3 
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which is 

($ 1; $ 2 ) = ($3,$!) = / ($ 3 ,$ 2 )=0. (; 
By differentiating these constraints and using ([!]), we have 

($1, ^ 1 ) x = ($ 2 , $ 2 )^ = ($3, $s)x = 

($!,$!), = ($ 2 , $ 2 ), = ($3,^3), = (: 

($1, $i) t = ($2, $2>t = ($3, $3)* = 
f x = i($ 3 , A$l) + i/(($ 3 , $3) " ($1, $l)) 

= i($ 3 , A$ 2 ) + i^(($s, $3) - ($2, $2)) 
= i($ 3 , A$l) + u# + i/(($ 3 , $3) - ($1, $1)) 
= -i($ 3 , A$ 2 > -uf- ig((®3, $3) - ($2, $2)) 



= 2i($ 2 ,A$ 1 ) +2i($ 2 ,^ 3 )($ 3 ,$ 1 ) 

" ($!,$!> "($2, $2) l " 

Mx = 2/<? (: 

Uy = , es ($i;$i) 2 _ 1 ($2;$2) [2i(^ A 2 $1 ) + ^($2, A$ 3 ) + 2i^($ 3 , A $1 ) 

+i^(2($ 3 , $3) - ($1, $l) - ($2, $2)) ( ' 

+u(M 2 - |/| 2 ) + u«$ 2 , A$ 2 ) - ($ 1; A$i))] 

Vl = 2|/| 2 V2 = 2|^| 2 . (1 
Therefore, the Lax pair (|l|) becomes 

$ M = iA$i + i/$ 3 $ 2 ,a- = iA$ 2 + ig$ 3 $3,z = + ig$2 ("- 

$ li2/ = iA$! + m$2 + i/$ 3 $2, y = - iA$ 2 - uj$ 3 

$ 3i2/ = i/$i - i#$ 2 (1 

$ 1)t = (-2iA 2 + i|u| 2 + m)^ + (-2uA + iu y )$ 2 + (-2i/A - 2/ y )$ 3 

$ 2 , t = (2mA + + (2iA 2 - i|w| 2 - iu 2 )$ 2 + (2i#A - 2g v )$ 3 (5 

$ 3 , t = (-2i/A + 2/,)$! + (2i^A + 2^)$ 2 - 2i(|/| 2 - \g\ 2 )$ 3 

where u, u y , t>i, i> 2 , /, g, f y , g y are given by (|I"2l) -(|l8r) respectively. 

Corresponding to the symplectic form (|IDD , the Poisson bracket of two functions £ and r] 
{(f) jk , i(p jk } is given by 

r> I = 1 V f ^ dT] d^__dr]_\ 

ii,Vi " 2 ^ ^Re(0 ifc ) dlm(</> ifc ) 01m(0 ifc ) <9Re(0 jfc )/ 

1 ^ / <9£ dr) <9£ dt] \ ^ 

i j£ \ d(f)jk d(f) jk d<j) jk dfak J ' 

Hereafter, we always look <fij k and <pj k as independent variables in differentiations. 



Theorem 1 The Hamiltonians for and $2l\ ) are given by 



H x 
H 1 



■($i,A$i) - ($ 2 , A$ 2 ) - |($i,$ 3 )| 2 - K$ 2 ,$ 3 )| 2 

-($1, A$i) + ($ 2 , A$ 2 ) - |($1,$ 3 )| 2 + |($2,$3>| 2 



(23) 
(24) 



+iw($i,$ 2 ) - m($2,$i) 
2($ 1; A 2 $!) - 2($ 2 , A 2 $ 2 ) + 4Re(($ 1 , $ 3 ) ($3, A$x)) 
-4Re(($ 2 , $ 3 >($ 3 , A$ 2 » + 2(($ 3 , $ 3 > - ($!, $i))($i, $ 3 )($ 3 , $i> 

-2(($ 3 , $3) - ($ 2 , $2>)($2, $ 3 >($3, $2) 

|($2, A$x) + ($ 2 , $ 3 >($3, $i)| 2 - 2Re(g($ 1; $ 2 )) (25) 



- ($ 2 , $2) 
respectively, which satisfy 

{H*, ($!, $ 2 )} = ($!, $ 2 )} = {#*, ($1, $2)} = 



(26) 



Here u is given by ( \Tb\ ) and q is given by 

The proof is obtained by direct computation. 

Therefore, we obtain Hamiltonian systems (^)-(25) with the constraint of Neumann type 
($1, $2) — 0. Any solution of the Hamiltonian equations 



l<ftjk,x — 
—i<Pjk,x 



dH x 

d<j)jk 
dH x 

d(f>jk 



i<Pik. 



- Wjk 



Jk,y 



dHy 



!) Jk 



ucp jk 
- i<f>jk,t 



(27) 



gives a solution of the DSI equation (g), where u, vi, t>2 are given by ([15f) and (J18|) respectively 



3 Integrability 

Now we consider the integrability of the Hamiltonian systems given by Theorem [I] on the 
submanifold 

S = {($ 1; $ 2 , $3, i$i, i$2, i$ 3 ) e lR 6iV I ($1, $ 2 ) = 0, ($ !, $i) ^ ($ 2 , $2)}. (28) 

Here we still use 3AT complex numbers and their complex conjugates ($1, $2, $ 3 , i5 ) 2? 1^3) to 
represent a point in IR 6Ar . Clearly, S has two connected components characterized by ($1, $1) > 
($2, $2) and ($2, $2) > respectively Since ^ (<& 2 ,$2), S. S is a 

QN — 2 dimensional real analytic manifold, on which the coordinates can be given by 3 i, • ■ ■ </> 3 tv 
(with their complex conjugates) and 2N — 1 of <f>u, • • • , 4>in-, 4>i\, " ' , 4>2N (with their complex 
conjugates) whenever the remaining one is non-zero. 
Define 

71 = Re($i, $ 2 ) = $2) + ($2, $i» 

\ (29) 

72 = Im($i, $ 2 ) = — (($1, - ($2, $1)) 
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then S is denned by two real- valued functions 71 = 72 = 0. Since 

{71,72} = - ($2, $2)) 



(30) 



is never zero on S, the symplectic form (|l0| ) on IR 6iV naturally induces a (nondegenerate) 
symplectic form on S. The corresponding Poisson bracket of two functions £, 77 on S 1 is still 
given by fl2|) if they satisfy {£, 7,} = 0, {77, 7i } = (j = 1, 2). 

Hereafter, the Poisson bracket { } always denotes the standard Poisson bracket fl22"|) on 
1R 6JV . 

Let 



L(A) 



V 



1 



N 1 
j=l A ~ X J 



■ 0:jO\, 02j01j O.yjOi j ^ 

01j02i 02j02j 0:i,0 2 ; 
V 02j03j 03j03j J 



(31) 



Then we have 



Lemma 1 L(X) satisfies the Lax equations 

L y = [V,L] L X = [U,L] L t = \W,L] 



(32) 



if and only if the constraints ( \Tq ) hold. 

Proof. Let J = diag(l, —1, 0), L = diag(l, 1, 0), <pj 



' then 



N 



L(X) = L + E 



x x 



Since V(Xj 



-V(X 



3)1 



N 



Ly{X) = E 



3=1 A - x j 

N -1 JV j 

E TTTTT^W'^i] " E ^^( A ) " ViXiUM 



3=1 A - A i 



i=i A ~ A J 

iV 



[7(A) J I(A)-Lo]-i[J ! EM*] 



Hence L„(A) = [K(A),L(A)] if and only if 



TV 



[l ,v(x)] = uJ2^] 

3=1 



(33) 



(34) 



(35) 



Written in the components, this is exactly the constraints ([12]). This proves that the first 
equation of (|32|) is equivalent to (0). When (|T2"D holds, the other two equations of (|32| ) are 
obtained similarly as the first one. The lemma is proved. 
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By Lemma [II trL k (k > 1) are all conserved. Expand trL fe as a Laurent series 



trL k = 2 + J2 



3=0 



(36) 



which is convergent absolutely and uniformly as |A| > maxi<j<N \Xj\, then all {Sj} are con- 



served. 



Moreover, we can show that any two of {£j } commute with each other. This follows from 
the following more general lemma. 

Lemma 2 Suppose M 2nr = {(q n , • ■ ■ , q ln , q 2 \, • • ■ , q 2n , • • ' , Qn, • • • , Qm,Pu, • • ' ,Pin, • • ' ,Pn, • • • , 
Pm)} ^ equipped with the standard symplectic form 



uJ 



dpj a A dq 



ja- 



(37) 



l<j<r 
Ka<n 



Denote qj = (qji, • • ■ , qj n ) T , Pj = (Pji, • • • ,Pjn) T - Let Ai, • • • , A„ be n distinct real numbers, A 
be an r x r constant matrix, 



( 



m(\)=a+y: 



i 



PlaQla P2 a qia ' " " PraQla 
Plaq2a £>2»<?2a ' ' ' PraQ2a 



\ 



(3* 



\ PlaQra P2a.qra ' ' ' Praqra / 

Then for any two complex numbers \, fi and any positive integers k, I, a, b with 1 < a, b < r, 



{tvM k (X),tvM l (fi)} = 

{trM fc (A), (p b ,q a )} = -k[A,M k - l (X)] ab . 



(39) 



Here (p b ,q a ) = T,a=lPbaQaa- 



Proof. The j-th row of ^ — is — (pi a , • ■ ■ ,p ra ) and the other rows are zero. Similarly, 

OQja X — X a 
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the i-th column of - - is — (qi a , ■ ■ ■ , q ra ) T all d the other columns are zero. Hence 

0p ja A - \ a 

i{trM*(A),trM'M} 



-,(M» ( ^) tr( M- W ^)- 

1 1 



E hrnr^*" 1 ^)^ -(m'- 1 ^))^ (40) 



a,b,j,a 



= l —Va a q ba [M k -\\),M l -\ii)} ab 

a,b,a A ~ Aa t* ~ Aa 

= E^-f— ^ VV-^fM^^A),]^'- 1 ^)]^ 

= -^tr(^M(\) - M^))[M k -\\),M l -\^ =0. 
This proves the first part. The second part is proved as follows. 



~{trM*(A),(p 6 ,g )} 

£ (v^M^A))^^ - T i r (M*- 1 (A)) a ,a p ta ') ( 41 ) 

l,a \ A ~~ A « A - A a / 

(M(A) - A)M fc_1 (A) - M k -\X)(M(X) - A)) 

= -[AtM^wu. 

The lemma is proved. 

From this lemma, we know that the set {£j k ^} is in involution. 

Remark 1 For the three-wave equation, [L&J wrote down such {£• }. It was proved that they 



lib 



3 

were in involution and 3N of them were independent. However, for the three-wave equation, 
the first term of L is diag(/3i, (3%, (3%) (f3\, fa, fiz are distinct) rather than diag(l, 1,0) here. In 
the case fl\ — 1, (3 2 — 1, 03 — 0, those {S^} are obviously still in involution. Hence we can 
also use the result of ]7^|/ to get the involution. But we prove it more directly and easily here. 
On the other hand, with the constraint ($ 1; $2) — 0, the independence of {£j } completely 
changes. Hence we should prove the independence in this constrained case. 
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For real A, suppose the eigenvalues of the Hermitian matrix L(X) are v%, V2 and z/3, then 
trL fc = T,j=i tf, while 



detQi-L(X)) = l[(n-u j ) 

= /i 3 - (Vi + V 2 + ^3)yU 2 + {V\V2 + + VsVl)^ - VxV2V z 

for any \i. Hence trL fc can be expressed by the coefficients of /1 in det(/^ — L(X)). 
Suppose L(A) = {L jk )x< jjk <z, then 



(42) 



V\ + V2 + Vz = trL 

+ ^2^3 + ^3^1 



^i^^ = det L 



Define 



trL 



00 c(0) 

£1 Xk+1 



E 



detL 



Lkj L 



kk 



(43) 



00 c*(2) 
\ -* c fc 



E 

l<i<fc<3 



Ljk 
Lkk 



00 

E fc 



(44) 



These } differ from the {-F^} in |TJ| with a multiple i. 

(i) 

To simplify the expressions of S k s, we define the following simpler but equivalent conserved 
integrals E^'s, which are nondegenerate linear combinations of £ k s. This also makes the 
expressions (|50|)-(|52D of the Hamiltonians and the proof of Theorem f| simpler. For m > 0, 
define 

E (l) _ o£(0) _ 



($ 1 ,A m $ 1 ) + ($ 2 ,A m $ 2 ) 



E E 

i<«<i<3 z=i 



(45) 



£■(2) , 
~ m+l "i" c m+l 



c(0) 



E 

z=o 



E 



($i,A'$ 2 > ($2,A m -'$ 2 ) i 

($1,A%) (^A^x) ($ 3 ,A fc $l) 

($i,A%) ($ 2 ,A^ 2 ) ($ 3 ,A fc $ 2 ) 
($1,^3) ($ 2 ,A^ 3 ) ($ 3 ,A fc $ 3 ) 

£7® =4 1) -4 0) = ($3,A-$3) 

(fci.A 1 " 1 ^) ($„A m -%> 



i+j+k=m—l 
i,j,k>0 



+ E E 

l<«<j<3 i=l 



(46) 



(47) 
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The above sums are zero if the upper bound is smaller than the lower bound. According to 
Lemma ^ we have 

Theorem 2 = and ($i,$ 2 )} = for any j,k = 1,2,3 and m,n > 0. 

Therefore, {E$J} are in involution on S. 

Define 

fi x = ($!, $j) ft 2 = ($ 2 , $ 2 ) fi 3 = ($3, $3>- (48) 

By (HD-®, 

= fix + fi 2 E (2) = n x Q 2 - !($!, $ 2 >| 2 E® = fi 3 - (49) 

Hence fli, fi 2 , fi 3 are expressed by E^'s and |($i, $ 2 )| 2 . 
The Hamiltonians in Theorem [I] can be expressed as 



H x 


= -eP - 




H y 


1 


+ 2Ef } 






l 


(2£«£« - 4E?> 


— O2 




+(4 1)2 - 


44 2) )^ + (E^ - 



i 1) E^-E^ 2 E^+2E^E^) 

(eW 2 -2EP)(eU + E®) 
f ) )(Q 1 -Q 2 )H v + (H x ) 2 -{H v ) 2 

2 1 A ( I / rfS rfS \ 1 2 1 I / /J\ rT\ \ I 2 \ I / *T» A\ \ I 2 



(50) 
(51) 



( M ($ 1 ,$ 2 )-M($ 2 ,$ 1 )) i + 4(K$ 1 ,$3)r + i($2,$3)r)i(*i,$2)r . (52) 



Each Hamiltonian is a function of ^ s and $ 2 )| 2 , ($1, < & 2 ) 2 , ($ 2 , $i) 2 , which is smooth 
near S. 

Let 71 and 7 2 be two real functions in (|2~9"D defining 5. For any two functions if i and H 2 on 
H 6JV with 7fc } = (j, jfe = 1, 2), define H x = H 2 if H x = H 2 and V#i = VH 2 on 5. In this 
case, Hi and H 2 have the same Hamiltonian vector field which is tangent to S. Hence they give 
the same Hamiltonian systems on S. Moreover, for any function F on JR 6N with {F, 7^} = 
(A; = 1,2), {H 1 ,F} = {H 2 , F} holds. 

From 



Q 1 -Q 2 = a-^E^ 2 - 4(4 2) + K$!, $ 2 )| 2 ) (53) 
where a = ±1 depends on the connected component of S 1 . Let 



K» = (Si, - S! 2 )ff» = -£<' ] E? + 2E<? - E^E™ - B< 1)2 B< 3) + 2£< 2 >E< 3) . 1 ' 

For 7! and 7 2 , only the quadratic terms appear in the expressions of H y and FL l . Hence we 
have 

H y = <L K y ( 55 ) 
H* = ^E^EP - 4^ 2) + 2(Ei 1)2 - 2E^)(E{ 1) + e[ 3) ) 

+ (Ej ) 1)2 - 4£< i 2) )# a; + (E^ - 2Ef ) )Ky + (^) 2 - ^2-(^) 2 ) • (56) 



10 



By (|50D , H x depends on E% s only. On the other hand, the right hand sides of (|55|) and 
(|56l) also depend on E^'s only. Therefore, we get the following result from Theorem [1] and 
Theorem 0. 

Theorem 3 Three Hamiltonians H x , H y , H l defined by Theorem [I] commute with each other: 

{H x , H y } = {H x , H 1 } = {H y , H 1 } = (57) 

and they satisfy 

{H x , ($!, $ 2 >} = ($!, $ 2 >} = {H\ ($!, $ 2 >} = 0. (58) 
Moreover, each E$ is conserved under the Hamiltonian flows given by H x , H y , H l respectively. 

Next we shall prove the integrability of these Hamiltonian systems. That is 

Theorem 4 3iV — 1 real-valued functions {Ejfi}} (j = 1,3; < m < N — 1) and {Effl} 
(0 < m < N — 2) are functionally independent in a dense open subset of S . 

Proof. Let 

51 = {($!, $ 2 , $ 3 , i$i, i$ 2 , i$ 3 ) e S | 0i )iV ^ 0, <p 1)N ^ 0, $i) > ($ 2 , $ 2 )} 

5 2 = {($1, $2, $3, i$i, i$2, i$ 3 ) e 5 I 01,7V ^ 0, 01,7V ^ 0, ($1, $1) < ($ 2 , $2)} (59) 
S = ^ U S 2 

then 5* is a dense open subset of S. Similar to S, S has also two connected components, which 
are Si and S2. In S, we can solve 02, /v, 4>2,n from the constraint ($1, $2) = as 

fo,N = ~ £ P-fci = ~ £ P-fai- (6°) 

J=1 017V J= i <PlN 

Hence Si has global coordinates 

6 = {0 y , i0i, (1 < j < N); 2j , i02, (1 < j < N - 1); 3j , i0 3 , (1 < j < iV)}. (61) 

Let P e Si be given by $1 = (1, 1, • • • , 1) T , $ 2 = e(l, 1, • • • , 1, -N+1) T , $ 3 = e(l, 1, • • • , if 
where e is a small real constant. Here $3 is chosen to be parallel with $1 so that the following 
computation will be simplified. Since 02/v and 02tv are functions of the variables in 0, we have, 
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at Pn 



Xf + 0(e) ^3^ = 0(e) 



d(j>ij 

_ m 

<90i,- 

•902,' 



0(e 2 



•902,' 



<9£« 

_ m 

•903,' 



0(6) 



0(e 



E 

1=0 



E 



77] 

^E 

«=0 



Ai + • • ■ + 

Ai + ■ ■ ■ + X N -i -(n- i)a<v x?~ l + {N- 1)XT~' 



'IN l=0 



($1,^) A^-^iiv 
($i,A<$ 2 ) AJTW 



+ 0(e 3 ) 
+ 0(e 3 ) 



(62) 



= iVe £ Ar'(Ai + • • • + X l N _, + Aj.) + 0(e 
dE^ J:l 9Pl 3 ) dM 3 ) 



0(i 



•902,' 



0(i 



•903,' 



Here the subscript j is taken from 1 to JV for <f>ij, 
checked that 



eX™ + 0(e d 



and from 1 to N — 1 for 



det E^(Ai + --- + Aj v _ 1 + Aj; 



\z=o 

Hence the Jacobian determinant 

,(1) 77(1) 



;)/ A' 1 IJ T?*- 1 -' 77(2) ; , • 1 - 1 1. ■'■->> /.•' 

u \ rj 1 " " " ; f^V-li -^0 ) • • • ) -^0 )" , )%-L 

5(011,- 
N N^N-1 



0<m<N-2 
l<j<N-l 



n(2) 77,(3) 



J = 



iV J] ^ - A*) 

l<j<k<N-l 



,(3) 



It can be 



(63) 



>1,N, H>2U ■ ■ ■ , <P2,N-1, (P31, ■ ■ ■ , (P3,N) 
2 



(64) 



n (Aj-a*) 

K l<j<k<N 



II (A,-A fc ) + 0( e 
i<i<fc<JV-i 



2N\ 



When e is small enough and e 7^ 0, J is non-zero near the point Pq. Since J is a rational 
function of 



{<f>lj, 101, (1 < j < AO! 02,, 102, (1 < J < AT - 1); 3j , 103, (1 < j < AT)} 



(65) 



J = identically if J is zero in an open subset of S±. Therefore, J is non-zero in a dense open 
subset of Si. Similarly, J is non-zero in a dense open subset of S 2 . Since all E^p are real- valued 

functions, the Jacobian matrix of (E^, • • - , pjv-i, -^o 2 ^, " ' ■> -^/v-2, Eo 3 \ • • • > ^jv-i) with respect 
to the real coordinates 



Re(0n), ■ • • , Re(0i >A r), Re(0 2 i), • ■ ■ , Re(0 2 ,jv_i), Re(0 3 i), • • • , Re(03 iiV ) 
lm(0n), • ■ ■ , Im(0i j7V ), Im(02i), • • • , lm(0 2 ,jv-i), Im(03i), • • - , Im(03,7v) 

is of full rank 3 A" — 1. The theorem is proved. 



(66) 



Theorem 5 The Hamiltonian systems given by Theorem [| are completely integrable on S in 
Liouville sense. 
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Proof. We have proved: (1) {E$} (j = 1, 2, 3; m = 0, 1, 2, • • •) are in involution on S (The- 
orem D. (2) {E$} (j = 1, 3; < m < N - 1) and {E$} (0 < m < N - 2) are functionally 
independent in a dense open subset of 5* (Theorem |4j). It remains to prove that the Hamiltonian 
vector fields of all {E^} are complete. This follows from the compactness of each level set, 
which is a closed subset of the compact set 

{($!, $ 2 , $ 3 , i$l, i$ 2 , i$ 3 ) G 5 | ($„ $,) = fy , j = 1, 2, 3} (67) 

where fljo (j = 1, 2, 3) are constants. Therefore, the Hamiltonian systems given by Theorem [l] 
are completely integrable 

4 Example: An explicit solution for N = 2 

A \ / 0j 

Now suppose iV = 2, A = , $ 7 = . Let R±, i?2, -R3, G and K be defined by 

^ 2 = l^f + l^| 2 (j = 1,2,3) 

G=|</> 1 | 2 + |0 2 | 2 + |0 3 | 2 -^ (68) 
K = (($i,A$i) + ($ 1 ,$ 3 )($3,$i))/ J R 2 

then from the above list of conserved integrals, we know that R\, R2, R3, G and K are all 
constants. Moreover, 

($ 2 ,A$ 2 ) + ($ 2 ,$ 3 >($ 3 ,$2> = R 2 2 (R 2 3 + \ + ix-K). (69) 

Let 

4>j = Rj cos 9j exp(ia,) = -Rj si n % ex P(ift) (70) 

then the constraint ($ 2 > ^i) — leads to 

6 2 = vQx + tt/2 + In exp(i(ft - ft - a 2 + «i)) = a (71) 

where / is an integer and a = ±1. Notice that ([7(]) is invariant under the transformation 
#1 — > — #i, ft — > ft + 7r. Hence we can always choose a = 1. 

Let 5 = ft — ft — a>i + «3 and p = cos 2 Q\. Substituting (|70| ) into the second equation of 
(|68|), we obtain 

R\ C os 2 9- A = G - (R\ - R 2 2 )p. (72) 

The third equation of ( |68"D gives 

A cos 2 6 1 ! + /i sin 2 6>i + Rj (cos 2 6»i cos 2 6» 3 + sin 2 Q x sin 2 6> 3 ) 
+2Rl cos 6*i sin 61 cos # 3 sin 6* 3 cos 5 = K 

from which 5 can be solved as a function of p. 
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(74) 



From the equation (|19D, we have 

9i,x — R3 cos $3 sin #3 sin 3 
#3,x = {Rl — R\) cos ^1 sm #i sm <5 

= A + i?| cos 2 6*3 + i?| cos # 3 sin 6*3 tan 8 X cos 5 
&2,x = A + i?g cos 2 #3 — R\ cos #3 sin 83 cot 6 1 ! cos 5 
®3,x = R\ cos 2 6 1 ! + R\ sin 2 #1 + (R\ — R%) cos #1 sin 81 tan #3 cos 5 
Pi,x = P + R3 sin 2 6*3 + R\ cos 83 sin 83 cot 8\ cos 5 
/?2,x = P + -R3 sin 2 6*3 — R\ cos #3 sin 6*3 tan 81 cos 5 
fc x = R\ sin 2 6 1 ! + R\ cos 2 6>i + (R 2 - i? 2 .) cos 6»i sin 8\ cot 6» 3 cos 5. 

The first equation of the above system leads to 

p x = — 2i? 2 cos 81 sin 6*i cos #3 sin 83 sin (75) 

Solving cos 5 from ([F3), we get the equation of p: 



Px = -^P(p) (76) 

where 

P(p) = 46(A - p)p 3 + (4(p - A)G + 4fc(p -if)-(A-p + & — a) 2 )p 2 



+ (4(K - p)G + 2(A - p + 6 - a) (AT + G - p - a))p - (K + G - p - a) 2 
which is a cubic polynomial, 



(77) 



a = Rj h = R\-R\. (78) 
Suppose b > 0, A > p and P has three different real roots pi < p 2 < P3. Moreover, suppose 

K +G- p-a^0 K -G - A+6^0 
max ( 0, ) < min ( 1, j 

Then the solution p can be expressed by elliptic functions of x. Let p = pi + (p 2 — Pi V 2 , then 



^ = ±p V /(l-^)(l_^2 ) (g0) 

where 

k = \j(p2 - Pi)/(ps - Pi) P = \A(A -p)(p 3 - Pi). (81) 
Hence u; = ± sn(p(x — £ )), 

p = pi + (p 2 - pi) sn 2 (p(x - z )) (82) 
where x is independent of x, but may depend on y and t. p is a periodic function of x. 
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Remark 2 Since pj is a root of P, ^FJj) leads to 



± Pj (l- Pj )(G-b Pj )(a-G + b Pj ) 

x2 (83) 
K — Xpj — p(l — pj) — pj(G — bpj) — (1 — Pj)(a — G + bpj) ) > 0. 

(This is actually equivalent to (fl^)). Hence, max(0, ^y 2 ) < Pi < P2 < min(l, ~) /io/ds i/ there 
is a solution locally, since < cos 2 81 < 1 and < cos 2 63 < 1 should be satisfied. This also 
guarantees that the solution is global because p\ < p < P 2- Moreover, under the assumption 
K + G-p-a^0andK-G-X + b^0, P(0) ^ 0, P(l) ^ 0. Hence < Pl < p 2 < 1 and 
< p < 1. 

Remark 3 Using the formulae 

1 - fc 2 sn 2 (0 = dn 2 (0 = ^ In 6(0 +C (84) 

where C is a certain constant, the previous solution p can be expressed by G function. 

In order to compute the y and t-equations, we first write down the expressions for u, f and 
g. They are 



2LR1P 

u 



exp(i(cti — ct 2 )) \{K — X — R\ cos 2 #3) cot Q\ 



I?2 p2 

— i? 2 cos # 3 sin #3 exp(i<5) 



exp(i(ai — a 2 )) i(p + R 3 sin 8 3 — K) tan #1 
2 V 

o cos #q sin 61 exn f-wn ( 85 ) 



+i? 3 cos #3 sin 6*3 exp(— i5) 
± P 2d2 

In 



4 /? / 

2 = 7W2 p2N2 ^3 + A + p - K) - Ap - AP 2 sin 2 fl 3 - cos 2 0< 



/ = -R1-R3 exp(i(ai — 0:3)) ( cos 6*3 cos#i + sin# 3 sin 6\ exp(io~)) 
g = R2R3 exp(i(a 2 — — cos9 3 sin^ + sin^ 3 cos#i exp(io~)). 

With the help of MAPLE, (p0|) and (^1|) are reduced to the following simple equations: 

01,2/ = 7l^l,x ^3,y = 7l^3,x 9l,t = J2@l,x ^3,t = 72#3,x 

2R\K 2R\(Rl + X + p — K) 

Oil,y = 7l«l,x - ^2 ^2 «2,y = 7l«2 j:E ^ ^2 

Iti — iX2 ix^ — 1X2 

ZRlRl ' (86) 

n x — n 2 

OCl,t = l2Ctl,x + C12 a 2) t = l20L2,x - C 2 1 «3,i = 72«3,z + C 3 

(Pj ~ a j)y = li(Pj ~ <*j)x (Pj ~ oij)t = j 2 (Pj - aj) x + C (j = 1, 2, 3) 
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where the constants 71, 72, Co, C\ 2 , C 2 \ and Co are given by 

r 2 + r\ 

2(Rj + 2 Rl - (R 2 + R 2 2 )(R 2 3 + A + p)) _ AR\R 2 2 {Rl + X + /2-2K) 
72 ~ Rl-BZ {R\ - Rlf 

2(\-v)(Ri + Ri) 



»7) 



f?2 p2 

^ = (Rj - i?|) 2 + 2jRj2jR ? ~~ ^ ~ V)( G L ^ — ~) 

2 D 2\/^ + / i ts\ \.. \nA , /id2 d 2j^/ d 2 



^0 - 02 r>2 

h — *U 

+ (R 2 - i? 2 )(^r^ - Ki) - Ap - Ai? 2 ) + 4RiR 2 K(Ri + A + n - K) 
-X(R 2 - R 2 )(Rt + R\) - 2R)K 2 \ = (1, 2) or (2, 1) 

K 1 = K K 2 = R\ + \ + n- K 

C3 = ^ r2 _ R2 y [rI(R\ + 3i?i - RlR 2 2 )(X + Rl) + r\(r\ - r\ + rIrDu 

-2{R 2 + Rl){R{ + R\)K - 2R 2 R 2 (R\ - R\ 



Hence 



P = Pi + (p2 ~ Pi) sn 2 (p(x + nf t y + 7 2 t - x )) (89) 

where xq is an arbitrary constant, p is given by (]8lf) and the parameter of the function sn is k 
given by (|8l|). 

The solutions of the DSI equation are 



iR\R 2 

u = 



x {{2K - A - fi - a + 6)p(0 + (p + a- fT-G)- iy/P(p(0)) (90) 
x exp y Q(p(0) + - 713/ + 72*) + KC12 + C 2 i)t 



and 



where 



t; 1 = 2i? 2 ((p-A)p(0 + ^-p) 
u 2 = 2i^((A-/i)p(0+if-^ + a) 

f = x + 712/ + 7 2 t - a; 

_ JZpf - iff (fl§ + A + - if) 
a ~~ i? 2 + R 2 

^ . . 26p 2 + (/I - 6 - 2G)p + (K + G - p - a) 

g(/?) = MW) 

/>(0 = Pi + (P2 - Pi) sn 2 (pO 

and the parameter k of the function sn is given by ([Sip . 



(91) 



(92) 
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u has no singularity when ([79| ) holds because in this case < p < 1. 

Suppose the minimal positive period of the function sn with parameter k is T(k) and 

A = W)l QM0)d( - (93) 

Then we have the following properties of the solutions: 

(1) u is a double periodic function on (x,y) plane. The period for x + jiy is T(k)/p, while 
the period for (A + a)x + (A — a)jiy is 2n. 

(2) u is periodic with respect to t if and only if — — . — - . is a rational 

(Ci 2 + C 2 i + A l2 + a 72 )T(A;) 

number. 

(3) |w| 2 , f i and t> 2 are periodic functions of x + 71?/ + 7 2 t only, and they extend constantly 
in a transversal direction on (x, y) plane. 

(4) The phase of u depends not only on the linear functions of x, y and t, but also on an sn 
function of x + jiy + 7 2 i. This can be obtained from ( |90|) and 

1 u x 2fg 

(&rgu) x = Re = Re 7^ constant (94) 

1 u m 

by using (|85|) and tedious computation. 

It is still interesting to solve more general periodic solutions using this method. 
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